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ABSTRACT 

We study the growth of galaxy masses, via gas accretion and galaxy mergers. We introduce a toy model that describes (in a single 
equation) how much baryonic mass is accreted and retained into galaxies as a function of halo mass and redshift. In our model, 
the evolution of the baryons differs from that of the dark matter because 1) gravitational shock heating and AGN jets suppress gas 
accretion mainly above a critical halo mass of M shoclc ~ 10 12 M G ; 2) the intergalactic medium after reionisation is too hot for accretion 
onto haloes with circular velocities u c j rc S! 40 km s" 1 ; 3) stellar feedback drives gas out of haloes, mainly those with v C ir C S! 120 km s -1 . 
We run our model on the merger trees of the haloes and subhaloes of a high-resolution dark matter cosmological simulation. The 
galaxy mass is taken as the maximum between the mass given by the toy model and the sum of the masses of its progenitors (reduced 
by tidal stripping). 

Designed to reproduce the present-day stellar mass function of galaxies, our model matches fairly well the evolution of the cosmic 
stellar density. It leads to the same z=0 relation between central galaxy stellar and halo mass as the one found by abundance matching 
and also as that previously measured at high mass on SDSS centrals. Our model also predicts a bimodal distribution (centrals and 
satellites) of stellar masses for given halo mass, in very good agreement with SDSS observations. 

The relative importance of mergers depends strongly on stellar mass (more than on halo mass). Massive galaxies with m slius > m cr j t ~ 
n b /fl m M shock ~ 10 U M G acquire most of their final mass through mergers (mostly major and gas-poor), as expected from our model's 
shutdown of gas accretion at high halo masses. However, although our mass resolution should see the effects of mergers down to 
Wstars — 10 10 - 6 /r'M G , we find that mergers are rare for ra stars < 10" h~ l M e . This is a consequence of the curvature of the stellar vs. 
halo mass relation set by the physical processes of our toy model and found with abundance matching. So gas accretion must be 
the dominant growth mechanism for intermediate and low mass galaxies, including dwarf ellipticals in clusters. The contribution of 
galaxy mergers terminating in haloes with mass Mhai < (thus presumably gas-rich) to the mass buildup of galaxies is small at 

all masses, but accounts for the bulk of the growth of ellipticals of intermediate mass (~ 10 I0 5 /7~'M Q ), which we predict must be the 
typical mass of ULIRGs. 

Key words. Galaxies: formation — Galaxies: evolution — Galaxies: luminosity function, mass function — Galaxies: interactions — 
Galaxies: ellipticals and lenticular, cD — Galaxies: haloes 

1. Introduction galaxies merging together mainly gas-rich (often called 'wet' 

mergers) or gas-poor ('dry' mergers)? How important are major 



. ^ ' In the standard theory d White & Ree si |1978fc |Blumenthal ej~iD mergers involving galaxies of similar mass, compared to more 
1 1 1984 IWhite & Frenkl |1991|), galaxy formation is a two-stage minor ones involving very different mass galaxies? What are the 
process. The gravitational instability of primordial density flue- respective r oles of feed back mechanisms such a s reionisation 



03 ; tuations forms dark matter haloes that merge hierarchically into Q f the IGM dReeslll986l) . supernovae explo sions (iDekel & Silkl 



larger and larger structures. Galaxies grow within these haloes [198(1 , jets from active galactic nuclei (AGN. ISilk & Reesl 1998b 
(i) by accreting gas that dissipates and falls to the centre of the anc j trie absence or presence of shock fronts between the in- 
gravitational potential wells of the dark matter and (ii) by merg- falling and virialised gas , depending on the mass of the halo 
ing with other galaxies after their haloes have become part of jBirnboim & Dekellflool ? How do the relative impacts of these 
larger structures. These two paths explain why there are two physical processes depend on the final (observable) stellar mass 
galaxy morphological types: in the standard theory, gas accre- f the galaxy? 

tion is the mechanism t hat builds up disc galaxies (i.e. spirals; To answer these questions, one must first estimate the 

|Fall&EfstafhiovJ|l980|), while galaxy merging is the mecha- build-up of dark matter haloes, e ither through Monte-Carlo 

rrism by which elliptical gala xies acquire the bulk of their mass merger-trees dLacev & Cold [l99l ISomerville & Kolattl \\99^k 



(|Toomre|ll977];|Mamon|[l99l. iNeistein & Deke]|l2008h or through cosmological N-body sim 



In the present article, we wish to revisit the question of how ulations, which also provide spatial information. If we assume 
galaxies acquired their present-day stellar mass. Was it mainly that the position of a galaxy is tracked by the centre of mass 
through mergers or through smooth gas accretion? Were the of its halo, then the halo merger rate will also give us the rate 
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at which galaxies merge (assuming the knowledge of the evo- 
lution of haloes once they become subhaloes of larger ones). It 
is more difficult to compute the stellar masses of galaxies prior 
to merging, which determine the contribution from gas accre- 
tion, because they depend not only on the gas mass that accretes 
onto galaxies but also on feedback processes that eject gas from 
galaxies. 



Coleetal. 1994 



There are two 
ity. Semi-analy tic 



approaches to deal 
model s (SAMs . 



with this complex 



also ICattaneo et al.l [2006 



Kauffmann et al 

— TT — 



1993; 



Croton et al.l 12006c ISomerville et alj 



Bower et al .112 006a: 
2001 lLo Faro et al.ll2009i 
and references therein) use simple recipes to describe the vari- 
ous physical processes affecting the baryons within haloes, and 
predict many observables. However, even the simplest SAMs 
quickly reach considerable complexity and involve a large num- 
ber of free parameters. 

The halo occupation distribution (HOD ) approach 
dBerlind & Weinberg! 120021: IConrov et al.l 120061: lYang et al l 
[2009F does not ma ke any assumption about the physics that 
pvern the numb er dBerlind & Weinberg] 20021). the l uminosity 
Yang et al.l2003l) or the stellar mass (lYang et al.l2009h of galax- 
ies within haloes. Instead, in the HOD app roach, one computes 
either the mean (abundanc e matching, see iMarinoni & Hudson! 
|2002[ IConrov et all 120061) or the statistical distributions that 
these properties must have as a function of halo mass and 
redshift to reproduce the observational data (assuming that the 
properties of dark matter haloes, i.e. the mass function and 
clustering, are modelled correctly). However, HOD models 
cannot tell us how much of this mass comes from gas accretion 
and how much comes from mergers. 

Our approach is intermediate between the two. We param- 
eterise, in a single equation, the mass of stars formed after 
gas accretion as a function of halo mass and redshift, m stars = 
Wstarst^haio, «ciic(A^haio, z)L where m stals takes into account vari- 
ous feedback effects that quench gas accretion and subsequent 
star formation. We compute galaxy mergers by following dark 
matter substructures. When the latter are no longer resolved, we 
merge the galaxies on a dynamical friction timescale using a 
formula calibrate d on cosmological hydrodynamic simulations 
dJiang et al.ll2008b . thus allowing the formation of galaxies with 

TWstars > m stars [M hllio , U c irc(Mialo, Z)]- 

Not only is our approach much simpler than fully fledged 
SAMs (we shall see that it contains only four parameters), it is 
also simp ler than "lighter" models such as the one recently pro- 
posed by iNeistein & Weinmannl (120101) . who parameterised the 
dependence of the time derivatives of the masses of the stellar, 
cold and hot gas components of galaxies as a function of halo 
mass and redshift. Admittedly, the simplicity of our approach 
limits its scope. But we believe that it makes robust predictions 
on the mass growth of galaxies, which provide a stepping stone 
to more complex analyses. At the same time, differently from 
HOD, our model allow us to follow mergers and thus to separate 
growth via gas accretion and growth via mergers. 

The structure of the article is as follows. In Section 2, we 
present the N-body simulation used to follow the gravitational 
evolution of the dark matter, and two simple formulae, one that 
expresses the timescale for orbital decay of unresolved sub- 
haloes, and the other that expresses the stellar mass at the centre 
of each halo if galaxies grew only by gas accretion. In Section 3, 
we present our results for the relative importance of gas accretion 
and mergers as a function of galaxy mass. Section 4 discusses 
our results and summarises our conclusions. 



2. The galaxy formation model 

2. 1. Gravitational evolution of the dark matter 

The hierarchical formation and clustering of dark matter haloes 
is followed by means of a cosmological A^-body simulation. 
The computational volume of the simulation presented here is 
a cube with side length L = 50/i _1 Mpc and periodic bound- 
ary conditions. The values of the cosmological parameters used 
to generate the initial conditions at redshift z - 50 by means 
of the Zel'dovich approximation dEfs tathiou et al.lfl985h and to 
run the simulation are those from the Wilkinson Microwave 
Anisotropy Probe's 5th-year analysis (WMAP5) combined with 
result s from type la supern ovae and baryonic acoustic oscilla- 
tions dKomatsuet al.ll2009h : Q m = 0.279, Q A = 0.721, Q b = 
0.046, h = 0.70, cr 8 = 0.817 and n s = 96. The simula tion was 
run with the tree-PM code GADGET2 dSpringelll20"05l) at a res- 
olution of 512 3 particles. We saved the outputs at 101 timesteps 
equally spaced in time (At = 92lr l Myr) from z — 10 to z — 0. 
The number of timesteps was chosen so that At is smaller than 
any merging timescale. The particle mass is 7.3 x 10 7 /z 1 M Q and 
the force softening is 2 h kpc. 

The MPI version of the AMIG/Q halo finder, AHF 
dKnollmann & Knebd 12009) PI was used to identify haloes and 
sub-haloes in each of th e saved outputsrj AHF is an improvement 
of the MHF halo finder dGill et al.ll2004IK which locates peaks in 
an adaptively smoothed density field as prospective halo cen- 
tres. For each of these density peaks, we determined the particle 
with the lowest gravitational potential and considered the sur- 
rounding particles sorted in distance to iteratively find the largest 
collection of bound particles encompassing the density peak, 
using their mean velocity as the kinetic reference. Only peaks 
with at least 20 bound particles are considered as haloes and re- 
tained for further analysis. Therefore, the minimum halo mass is 
1 .5 x lO 9 /!^ 1 M Q (corresponding to D c ; rc = 15 km s~'). Our halo 
finding algorithm automatic ally identifies haloes, sub-h aloes, 
sub-sub-haloes, and so on (see lKnollmann & Knebd2009l for the 
details of the algorithm). 

For each halo at epoch z, we compute the virial radius r v ; r , 
which is the radius where the mean density drops below A(z) 
times the critical density of the Universe at redshift z. The 
threshold A(z) is computed using the spherical top-hat collapse 
model and is a fu nction of both cosmological model and time 
(iNakamurall 19961 cited in iKitavama & Sutolll996t lGrosslll997l: 
iBrvan & Normanlll998l) . For the cosmology that we are using, 
A = 98 at z — 0, i.e. the overdensity at the virial radius is 352 
times the present-day mean density of the Universe. 

After each halo has been identified by AHF in this manner, 
we define sub-haloes to be haloes (defined with the same over- 
density A) that lie within the virial region of a more massive 
halo, the so-called host halo, of mass M^Sl- As sub-haloes are 
embedded within their host halo, their own density profile usu- 
ally shows a characteristic upturn at the radius r t < r v j r where the 
density profile becomes dominated by the surrounding host halo 
FJ. We use this truncation radius as the outer edge of the sub-halo. 
Halo and sub-halo properties (i.e. mass, density profile, velocity 



Adaptive Mesh Investigations of Galaxy Assembly 

2 AMIGA'S Halo Finder 

3 AHF can be downloaded from http : //popia. ft .uam. es/AMIGA 

4 The actual density profile of sub-haloes after the removal of 
the host's background drops faster than for isolated haloes (e.g. 
iKazantzidis et alj|2004l) . Only when measured within the background 
still present shall we find the characteristic upturn used here to define 
the truncation radius r t . 
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dispersion, rotation curve) are calculated using the gravitation- 
ally bound particles inside either the virial radius r vn - for a host 
halo or the truncation radius r, for a sub-halo. 

We build merger trees by cross-correlating haloes in consec- 
utive simulation outputs. For this purpose, we use a tool that 
comes with the AHF package, called MergerTree. As the name 
suggests, it serves the purpose of identifying corresponding ob- 
jects in a same simulation at different redshifts. We follow each 
halo (either host or sub-halo) identified at redshift z-Q back- 
wards in time. The direct progenitor at the previous redshift is 
the one that shares the greatest number of particles with the 
present halo and is closest to it in mass. The latter criterion is 
important for sub-haloes as all their particles are also part of the 
host halo, but there is normally a large gap between sub-halo and 
host halo masses. 

In the study presented here, we walk along the tree start- 
ing at redshift z = 10 and moving toward z-Q. Sub-haloes are 
followed within the environment of their respective hosts until 
the point where AHF can no longer resolve them either because 
they have been tidally disrupted or because they have merged 
with their host (of course, many sub-haloes survive until z=0). 
Section|22]below describes how we follow sub-haloes when we 
can no longer resolve them numerically. 

2.2. Growth via galaxy mergers 

We associate galaxies to the smallest resolved dark matter sub- 
structures. We say that two haloes have merged at the cosmic 
time t^ loes when a subhalo is no longer resolved in the N-body 
simulation (Fig. 1). With this definition, the merging of two dark 
matter haloes is a necessary but non-sufficient condition for the 
merging of their central galaxies because generally a sub-halo 
ceases to be resolved before reaching the centre of its host halo. 
Let fdf be the time the sub-halo takes to reach the centre after it is 
no longer resolved. The galaxy at the centre of the sub-halo has 
merged with the galaxy at the centre of the host halo by cosmic 
time f if fj^ 10 " + fdf ^ t. If this condition is satisfied, then the 
subhalo's central galaxy is merged with that of the host halo at 
the first timestep for which t > fJ}f loes + fdf - Otherwise, the galaxy 
properties are frozen after f^ aloes . 

To compute the time for the satell ite's orbit to decay to the 
centre, we use Chandrasekhar's ( 1 1943b dynamical friction time, 
which can be written as 



f df = A 



(1) 

GM sat lnA' v ; 

(e.g. lBinnev & Tremaind20"08l) . where r is the sub-halo's orbital 
radius at the time t^ loes when it ceases to be resolved, i> c ; lc is 
the circular velocity at radius r, InA is the so-called Coulomb 
logarithm, a nd A is a factor of order unity that d epends on orbital 
eccentricity dColpi. Maver. & Governatolll999L and see below). 
In the following two paragraphs, we explain our choice for In A 
and A. 

Ijiang et all (120081) used high-resolution cosmological hydro- 
dynamical simulations to show that, once a halo of mass M sat 
penetrates the virial sphere of a larger one, of mass Mh a i (and 
thus becomes a sub-halo, i.e. satellite), the time it takes to fall 
to the centre of its parent halo is given by equation ([TJ, where r 
and u c j rc are measured at the virial radius r v j, of the host halo and 
where the Coulomb logarithm is In A = ln(l + M^ a \ /M s . dt ). We 
thus assume that for r < r v ;,-, the appropriate Coulomb logarithm 
is the one computed using 

A=l + ^2, (2) 

M sa , 




Fig. 1. Sketch illustrating the merger of two haloes. Once the 
smaller halo penetrates the larger one (whose central galaxy is 
denoted by a red circle), it bec omes a sub-halo. It s orbit (adapted 
from an A^-body simulation bv lColpi et al.l (119991) ). shown by the 
blue curve (with positions shown for equal timesteps), decays by 
dynamical friction and its size, shown by the circles, is reduced 
by tidal stripping from the larger halo's potential. The solid part 
of the blue curve is followed directly by the N-body simulation, 
while the dashed part shows the trajectory that the sub-halo is 
expected to take after it {open circles) is no longer resolved in the 
simulation. In our analysis, the duration of the motion along the 
dashed stretch of the blue curve is computed analytically from 
Eqs.tQ} and ©. 



where Mh a \ (r) is the halo mass enclosed within the sphere of 
radius r. 

The normalisation A of the dynamical friction time is often 
chosen to be inde pendent of the shape of the satellite's orbit (e.g. 
1.17 according to lBinnev & Trema ine 2008). But for given orbit 
apocentre, the time of orbital decay is expected to be shorter for 
satellites on radial orbits, which feel more dynamical friction as 
they reach the high density region near the centre of their halo. 
iJiang et all calibrated this effect and found 



A(e) = 1.17 (0.94e 6 + 0.6), 



(3) 



where e = J/J CUC (E) is the orbital circularity, that is, the ra- 
tio of angular momentum to that of a circular orbit of the same 
ene rgy (e - for a radial orbit and e = 1 for a circular or- 
bit). |jiarig_eTal]'s parametrisation of the Coulomb logarithm and 
normalisation of the orbital decay time, based upon hdyrody- 
namical cosmological N-body simulations, implicitly takes into 
account the greater concentration of baryons with respect to the 
dark matter as well as the effects of tidal stripping on the dynam- 
ical friction time. 

We now explain how we compute in our model the value 
of e that we insert into Eq. (f3]). The radial coordinate r, the 
orbital speed u, and the radial and tangential velocities v r and 
v, — yjv 2 - v 2 of the sub-halo with respect to the centre of mass 
of the host halo are extracted from the N-body simulation at the 
last timestep when the sub-halo is resolved. We compute the sub- 
halo's angular momentum and energy per unit mass. They are, 
respectively, J — rv t and E = 1/2 u 2 + 4>{r), where (p(r) is the 
gravitational potenti al of the host halo th at we compute assum- 
ing an NFW model dNavarro et al 1 ll996h . This requires to know 
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the halo concentration. The latter could be evaluated directly 
from the particle distribution but this may lead to noise, particu- 
larly in low-mass haloes, which contain fewer particles. Instead, 
we use the most recent determination of the concentration- 
mass relation for regular haloes in ACD M simulations wit h 
WMAP-5th year cosmological parameters dMaccio et alj|2008b : 
c - 6.1 (M ha i o /lO 14 /iM ) -0 094 . With this assumption we com- 
pute J C uc(E) by solving the equation E = l/2i?- + <p{r) for 
r where the squared circular velocity is y?. = GM^°?Hr)/r, 
and where M?°j St (r) is also computed with the NFW model. This 
equation has a solution for E < 0. When E > 0, the sub-structure 
is not gravitationally bound to the host, so the satellite and the 
central galaxy do not merge. 

In our analysis, the galaxy merging history is determined 
fully and solely by the underlying evolution of the dark mat- 
ter component. It does not depend on the efficiency with which 
gas dissipates, sinks to the centre and makes stars. This is why 
we have discussed galaxy mergers immediately after presenting 
the N-body simulation. The processes that determine the growth 
of luminous galaxies via accretion within dark matter haloes are 
considered in the following subsection. 



2.3. Growth via gas accretion 

We compute the stellar mass of a galaxy, m stars , as a function of 
Mhaio usi ng a model that is an improvement over the one intro- 
duced by ICattaneol (l200 1 1) . Despite its simplicity, that model was 
shown to be in reasonable agreement with the galaxy luminos- 
ity function and the cosmic evolution of the star formation rate 
(SFR) density. 

The key assumption is that, in the absence of mergers, the 
mass of the stars in a dark matter halo is basically a function of 
two quantities only: the halo mass Mh a i and the halo circular 
velocity y circ (M ha i ,z). 

The effects of shock heating and feedback enter our model 
through the form of this function. Cosmic reionisation sup- 
presses gas accretion and star formation in haloes below a min- 
imum circular velocity u re j on , stellar feedback mitigates star for- 
mation in haloes below a characteristic circular velocity usn, an d 
gravitational shock heating coupled to black hole feedback sup- 
presses gas accretion and star formation above a critical halo 
mass of M s hock- 

The stellar mass of a galaxy at any given time is assumed to 
be given by 



Mstars = max jm stars [Mh a i ,y c irc(M halo ,z)], ^ 

{ prog 



'"stars ^prog | 



(4) 



where 2 prog * s me sum °f me stellar masses of all pro- 

genitors and < Sp r0 g < 1 is a factor that accounts for the tidal 
stripping of the mass of a satellite progenitor prior to merging. 
For the main progenitor, s plog = 1. The way in which s pmg is 
computed for the other progenitors is discussed in Section l2%l 

Our model does not contain any rate and does not integrate 
any differential equations. However, Eq. (|4]i implies that the stel- 
lar mass at timestep ; can be larger than the value computed from 
the function m stals at timestep i if the stellar mass at timestep 
r — 1 is larger than the mass computed from the function m stals 
at timestep i. This has a simple physical reason: besides the ef- 
fects of stellar evolution (neglected in our model), the mass of 
the stars that have already formed and are already in a galaxy 
cannot decrease, even if the potential well becomes shallower. 



While equation suggests that the final halo mass is in- 
dependent of the mass accretion history (MAH) of its halo, our 
prevention of any decrease in stellar mass implies some depen- 
dence on the halo MAH: A galaxy with a low-mass present-day 
halo might have never formed stars if its halo circular velocity 
was always below the threshold for galaxy (star) formation (see 
Sect. 12.3.21 below). Another galaxy with the same halo mass at 
z — may have acquired stars if its halo mass had grown fast 
enough for its circular velocity to break the threshold for galaxy 
formation, and in our model it would keep its stellar mass. The 
implications for the results of our model will be discussed in 
Section 3. 

We now describe how we choose the functional form of ra stars 
(the rest of this Section) and s pmg (Section l24t . 

2.3.1. Shock heating and AGN feedback 

The accretion of gas onto low mass h aloes is understoo d to 
proceed in the form of cold filaments dKeres et al.l F2005) that 
fall onto the disc, where the new material is shock heated, but 
cools very rapidly thanks to the very high gas density of the 
disc. However, as shown by iBirnboim & Dekel d2003l) . above 
a critical halo mas s M^ocM ~ \0 n M n (roug hly independent in 
time, see Fig. 2 of iDekel & Birnboimll2006l) . the condition for 
the propagation of a stable shock is satisfied, and the post-shock 
gas builds up a h ot atmosphere that cools inefficiently (see also 
lOser et al Moreover, AGN couple to the h ot gas prevent- 

ing it from cooling down again (see the re view of Cattaneo et al.l 
l2009l and references therein, in particular lCroton et al.|l2006 | for 
the firs t explicit implementation in a SAM). Following ICattaneol 
(l200lh . we model these two effects by assuming that the mass 
accreted by the galaxy from the start of the simulation to the 
timestep under consideration is, to first order, 



1 + M Mo /M sh ock 



(5) 



for v clIC » Ureion, where f D = £2b/£^m is the mean baryonic 
fraction of the Universe. In other words, m a J? cr - ft, Mhaio when 
M ha io « M shock , while 

"^accr fb M S hook when Mhdo » M shock . 
In the latter case, the galaxy mass grows very slowly with in- 
creasing halo mass (except for the effects of mergers, which are 
not considered in the calculation of the accreted mass)0 



2.3.2. Reionisation 

Heating by the photoionising UV background raises the 
entropy of the gas and suppresses the con c entrat i on of 
baryo ns in shallow potential wells ( Ikeuchil 19861: iReesI 



1986; 



iBlanchard. Valls-Gabaud. & Mamon 
HEfstafhiouT 



Compton heating - , 
iKravtsov, Gnedin. & Klvpinl |2004|) 



1992; 



19921 — for 



Here 



Gnedinl 12005 
take a phe- 



nomenological approach and model the effects of photoioni- 
sation heating by introducing a lower circular velocity cut-off 
at Voire ~ ^r eion that mimics the shutoff of gas infall whe n 
r v ir < r IGM (IBlanchard et al.lll992t Ffhoul & Weinberg! 1 19961) . 
Cosmological hydrodynamical simulations give conflicting 
answers on the temperature of the intergalactic medium (IGM) 
in the intermediate density regions outside the vi rial radius, from 
which gas should fall in JMcO uinn et al 1 I2Q09L and references 



5 Admittedly, the denominator of equation (T5J is ad hoc. We tried 
powers other than unity and other forms, but the form of equation lf5J 
was the simplest one that gave us a good match to the high-end of the 
mass function (see Sect. 13. It . 
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therein). We assume that this IGM temperature does not vary 
with redshift at z ~ z re ion > 10 (z — 10 is where our calculations 
start), whic h is in rough agree ment with hydrodynamical 
simulations dMcOuinn et al.l 12009:, and references therein). We 
also suppose that the mass m accr that can flow to the centre in a 
halo of mass Mh a i is suppressed with respect to the first-order 

accreted maSS by a factor Of 1 - (v re i on /v c i rc ) 2 at U c i rc ^ Zreion, so 

that: 



(6) 



At Ucirc < freion the suppression is total (m accr = 0). This 
extreme approximation is ad hoc, but it is not unreasonable 
given the high mass-to-light ratios measured in low-mass ob- 
j ects and is in agreeme nt with hydrodynamical simulations 
(iThoul & Weinbergllio-o-fih l 5 ! 



2.3.3. Supernovae 

When the circular velocity becomes larger than u re j on , the gas 
starts flowing to the centre and making stars, which feed en- 
ergy back to the interstellar medium mainly via type II super- 
nova explosions. The energy fed back to the interstellar medium 
is proportional to the stellar mass of the galaxy, ra stars . A fraction 
of this energy is used to drive an outflow. Let wi w ind be the total 
mass ejected from a galaxy from the start of the simulation to the 
timestep under consideration. For a given wind kinetic energy, 
1 /2 m w in,iy 2 ind , OT w i n d is maximum when the wind speed, u w ind, is 
minimum. However, the wind speed cannot be lower than the es- 
cape velocity, which is proportional to the circular velocity u c i rc , 
because otherwise the gas would not flow out. The maximum 
outflow condition t hus implies that m w ; n d = (vs^ / i! c i rc ) 2 m stm - s 
(Dekel & Silk 1986), where vl N is a proportionality constant, the 
value of which is to be determined by fitting the galaxy mass 
function and is physically related to the supernova energy that 
goes into wind kinetic energy. We also require mass conserva- 
tion: OT w ; n d + m stars = m accr . In other words, all the gas accreted 
onto a galaxy is either ejected or forms stars. This assumption 
is crude because galaxies do not turn all the accreted gas in- 
stantaneously into stars. However, even in spiral galaxies, typi- 
cal gas fractions (cold gas divided by cold gas plus stars) rarely 
exceed ~ 10 - 20% (except for dwarf irre gulars and blue com- 
pact dwarfs, see e.g. lMcGaugh et al1l2010l) . Therefore, the error 
on the galaxy stellar mass that one makes by assuming that the 
entire galaxy mass is in stars is consistent with the other un- 
certainties of our model at all but the lowest masses. When we 
substitute m wind = (usN/fcirc) 2 Mstars into ra wind + m stars = OT accr , 
we find 



(7) 



2.3.4. Synthesis 

By putting together Eqs. (0, (O, and ||7), we can write a simple 
but physically motivated equation for the model stellar mass (in 



6 iGnedinl d200Ch have used more refined hydrodynamical simula- 
tions to find that the suppression by reionisation scales as u' irt . (see also 



Okamoto et al. 2008), but this should produce negligible differences in 
the mass range of galaxies studied here. 



the absence of mergers) of the central galaxy of a halo of mass 
Mhzio in the absence of merging: 



fbM h 



alo 



1 + M halo /M shock 



(8) 



where 



=(M h alo,Z) 



A(z) 



1/6 



[GH{z)] lli M, 



1/3 jLfl/3 



halo 



(9) 



Had we taken into account that some gas is neither blown away 
by the supernovae nor converted into stars, then mass conserva- 
tion would have given m w j nd + m stars /(l - g) — m accr , where g 
is the fraction of galaxy baryons in the form of remaining gas. 
Hence, equation (0 would become 



fb A/halo 



1 + Mhaio/Mshock 



(10) 



At each timestep, the galaxy stellar mass is assigned accord- 
ing to equations (0 and (0, unless this produces a lower stellar 
mass than at the previous timestep. 

Equation ([8]) should be seen more as an empirical fitting for- 
mula than as a physical model, although the arguments that we 
have presented above provide some physical justification for it. 
We are neglecting a host of other physical processes such as the 
conditions for the formation of clumpy molecular clouds from 
cold atomic gas, and mechanisms operating inside halos such as 
ram pressure stripping, conduction, magnetic fields, cosmic rays, 
etc., because we wish to build the simplest model that provides 
an adequate fit to the galaxy mass function of galaxies and thus 
allows the study of the mass buildup of galaxies. 

The key point of this article is not to present a new model 
for processes such as shock heating, cooling, star formation, and 
feedback (SAMs do that more satisfactorily), but rather to esti- 
mate the role of mergers in the galaxy growth. For this we need 
a prescription to associate a galaxy mass to a halo mass at each 
redshift and this is what Equation (0 does. What is more im- 
portant is that we compute galaxy mergers accurately and this is 
where our model is state-of-the-art (Section l2~2l i. 

Figure [2] illustrates how, in our toy model (eq. (8j), the ef- 
ficiency with which halo mass growth results in stellar mass 
growth depends on halo mass. This efficiency has been plot- 
ted, at two different redshifts, for the parameter set that gives 
the best fit to the z — galaxy mass function (see Sect. 13.11 and 
Table Q] below). The star formation efficiency is maximum at 
Mhaio ~ 10 11 - 10 12 /2~' M Q , decreasing sharply at lower masses 
and more mildly at higher masses. Figure [2] shows that in a z=0 
halo, < 20% of the baryonic mass is expected to be in stars. The 
fraction of baryons that may be locked up in stars can reach 50% 
at high redshift, where feedback is less important because the 
circular velocity thresholds of equation (O correspond to lower 
mass thresholds (eq. Q, noting that A(z) and especially H(z) 
increase with z). Our model produces a z = stella r mass vs. 
halo m ass relation very similar to that obtained by iGuo et all 
d2010h with the abundance matching technique, except that our 
toy model (eq. O) underestimates the stellar mass at large halo 
masses relative to the abundance matching solution. This will be 
compensated by the effects of galaxy mergers (see eq. @) as we 
shall see in Sect. [3](Fig. [6]) below. 

Figure [3] shows how these processes quench star formation 
during the growth of a dark matter halo. The contours show 



5 



A. Cattaneo, G. A. Mamon, K. Warnick, A. Knebe: How do galaxies acquire their mass? 



log M halo (h" 1 if.) 
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50 100 500 1000 
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Fig. 2. Star formation efficiency versus halo circular velocity and 
mass at two different epochs in our toy model (eq. O). This 
figure has been plotted for the parameters of Table[TJthat produce 
the best fit to the stellar mass function at z — (see Sect. 13.11 ). 
The upper (black) mass scale is for z=0, while the lower (red) 
mass scale is for z=5. The gold shad ed region shows th e relation 
predicted by abundance matching bv lGuo et alj d2010h . 



m stai-s/(/b A^haio) as a function of halo mass and redshift (for the 
best-fit parameters used for Fig, [2] see Sect. 13.11 and Table [TJ. 
The term in y c ; rc in Eq. © suppresses galaxy formation in low- 
mass haloes, while the cut-off at M s h oc k suppresses galaxy for- 
mation in high-mass haloes. The result is a galaxy formation 
zone at 10 10 M o < M hnio h £ 10 I25 M o . Galaxy formation begins 
and ends when a halo moves in and out of the galaxy forma- 
tion zone, respectively. The four red curves in Figure|3]show the 
mass growth with redshift of four haloes whose final masses at 
z=0 are roughly 10 11 , 10 12 , 10 13 and lO 14 /2-'M . The haloes that 
are more massive enter and exit the galaxy formation zone at an 
earlier epoch. Since the most massive haloes contain the most 
massive galaxies, this explains why the most massive galaxies 
con tain the oldest stellar populations (archaeological downsiz- 
ing, ^|CattaneoietaD|2008|). 

The dashed curve in Figure [3] shows the characteristic mass 
Mni of density fluctuations that become non-linear at redshift z, 
i.e. the characteristic mass scale on which the cosmic variance 
of linearly extrapolated primordial density fluctuations is equal 
to (5 C - 1 .68 at redshift z. Here 5 Q is the density contrast of a top 
hat fluctuation that the linear theory predicts at the time when 
the fluctuation actually collapses. The value of 8 Q depends very 
weakly on O m and Qa. The redshifts at which M n \ comes in and 
out of the galaxy formation zone mark the beginning and the end 
of the main epoch of galaxy formation in the Universe. 

Similar simple toy models were independentl y constructed 
by D. Croton (private communication) and by iBouche et al.l 
(1201 Ol) . although in their models the low-mass cut-off does not 
depend on redshift, while in our model it is the circular veloc- 
ity that is fixed, so the minimum mass decreases with increasing 
redshift (see Figure|3]l. 
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Fig. 3. Growth of halo mass and efficiency of star formation. The 
contours show the efficiency of star formation, m^^/ifi, MhaioX 
as a function of redshift and halo mass, for the set of parame- 
ters that fit best the galaxy mass function (Sect. |3~D and Table[TJ. 
The horizontal dash-dotted line shows the critical halo mass for 
shock heating, Mh a i = Mshock- Baryons make stars efficiently 
only in the horizontal band of the diagram shown by the shaded 
region. The red irregular lines show the mass growth of four 
haloes, whose present masses are about 10 11 , 10 12 , 10° and 
10 14 /z _1 M o . They are backwards tracks for the most massive 
progenitors. The diagonal dashed line shows the characteristic 
mass scale M n \ on which density fluctuations become non-linear 
at redshift z, i.e. where cr(M n \,z) = 1.68. The horizontal green 
line shows the mass resolution M res = 1.5 x 1O 9 /i~'M of the 
N-body simulation. 



2.4. Tidal stripping 

The dynamical fr iction time given b y eqs. ([TJ, (f2]i and 0), which 
was calibrated bv lJiang et al.l (l2008h on hydrodynamical cosmo- 
logical simulations, implicitly incorporates the effects of tides, 
which strip the secondaries to lower masses as they orbit through 
the primary, on the dynamical friction time. However, knowl- 
edge of the stellar mass at a given time is necessary to compare 
to the observed galaxy mass function, so we need to estimate the 
time evolution of the stellar mass caused by tidal stripping. 

Tidal forces are known to strip galaxy hal oes very ef- 
ficien tly once they penetrate larger ones (e.g. IGhigna et al.l 
Il998[). and each per icentric passage generates more mass loss 
(lHavashi et aT1l2003h . But tides also affect, t o a lesser extent, the 
more bound stellar material. For example, iKlimentowski et all 
(I2009h ran simulations of a low-mass high-resolution spiral 
galaxy around the fixed potential of the Milky Way and found 
that over 5 orbits, while the dark matter mass was reduced by a 
factor 100, the stellar mass was reduced by a factor 10, meaning 
that a fraction of 77 str ; p = 1 - (1/1 0) 1 ^ 5 = 0.37 of the stellar mass 
was lost at every pericentric passage. This tidally stripped stellar 
mass should form what is known as the stellar halo in the case 
of a galaxy and the intracluster light in the case of a cluster. 

According to the simulations of IKlimentowski et all the 
mean mass loss 77 stl i p is roughly constant for every pericentric 
passage, so the importance of tidal stripping depends on the 
number of orbital revolutions that a galaxy makes from when 
its orbit starts decaying to when either the galaxy merges with 
the central galaxy of its parent halo or the simulation ends. In 
the former case, the total time span during which a galaxy is 
stripped is given by the time it takes to decay by dynamical fric- 
tion from the virial radius to the centre (Eq. [TJ with r and t; c i rc 
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Table 1. best-fit parameters corresponding to the full model in 
Figure [4] 



^reion 


40kms -1 




120 kms" 1 


M shock 


8 x W n hT 1 M G 


77strip 


0.4 



taken at the virial radius). Taking f or b 
time, we find 



2nr/v C a C for the orbital 



T = 



^orb 



A(e) ^vir 

~2n GM si 



In A 



Me) 

2n ln(l + MS/M sat ) 



(11) 



The normalisation, A(e), o f the dynamical fr iction time is com- 
puted using the formula in lJiang et al.l (120081) . which fits the re- 
sults of SPH simulations and accounts for the tidal stripping of 
sub-structures. Eq. (fTTT i implies that, for a fixed value of ?7 str ip, 
the importance of tidal stripping increases with M^^/M SRt . 



2.5. Summary 

In summary, the final equation for the galaxy stellar at any given 
time is 



max < 



"'stars 



i 

Piog 



'"stars v 1 



?7strip) 



(12) 



where m st ars and r are given in equations © and (TTTI) . respec- 
tively, with t — for the main progenitor. Equation (fT2l con- 
tains four free parameters (u re i on , usn, Mhock, and ^strip), the first 
three of which enter the model through the function m stars of 
equations dS) and ©. The fourth parameter 77 stl ; p affects the fi- 
nal galaxy masses by controlling the efficiency of tidal stripping, 
but, as we shall see, n str i p has comparatively little bearing on our 
results. There are no free parameters in the calculation of the 
halo and thus galaxy merger rate. 



3. Results 

3. 1 . Best-fit model and effects of physical processes 

We start with the set of our 4 parameters, v le i on , v$n, Mshock, and 
77 st ri P , that p rovide a b e st-fit t o the galaxy stellar mass function 
deduced bv lBell et all (l2003h from the SDSS observations, and 
then consider how the predictions vary when the assumptions are 
changed. 

The stellar masses of galaxies derived by iBell et all d2003l) 
to compute the SDSS mass function depend on the Hubble con- 
stant used to compute galaxy luminosities and on the stellar ini- 
tial mass function (IMF). The masses inferred from observed 
luminosities scale as lr 2 while simulated masses scale as h~ x . 
As it is not possible to make the comparison between model and 
observations completely independent of h, we convert the ob- 
servatio nally determined masses in units of h~ l M Q by assuming 
h = 0.7. lBell et al.l d2003h compute the stellar m ass function us- 
ing a diet Salpe ter I MF (IBell & de JondfeOOlh . To convert to a 
ISalpeterl d 19551) or a lKroupald2002h IMF one should add a cor - 
rection of +0. 15 and -0. 1 dex, respectively (see lBell et al.l2003l) . 

The solid lines in the two panels of Figure [4]indicate th at our 
best-fit model (Table Q} matches well the lBell et alJ (12003) mass 
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Fig. 4. Stellar mass function of galaxies at z-0. The filled 
squares show the galaxy mass function determined observation- 
ally by Bell et al. (2003). The solid line shows our best-fit model 
(the 'full' model, with parameter values listed in Table [TJ. The 
fit assumes h = 0.7 because observationally inferred masses 
scale as hr 1 . (a) Variations of the stellar mass function for dif- 
ferent assumptions on merging and stripping: best-fit full model 
(black solid line), no tidal stripping (?7strip — 0, red dotted line), 
no mergers at all (green dashed line), and maximum merging 
(blue dash-dotted line), (b) Variations of the stellar mass func- 
tion for different assumptions on the baryonic physics: best-fit 
full model (black solid line, same as in panel a), baryons trace 
dark matter (m stars = /b^haio, red dotted line), no high-mass 
cutoff (M s hock — » °°> green dashed line), and no low circular ve- 
locity cutoff (i>sn = freion = 0, blue dash-dotted line). In both 
panels, all variations start from the best-fit parameters of the full 
model except for the fact that we use 77 stl j p = in all models apart 
from the full one. 



function at 10 M a < hm st . dIS < 10 M Q . Differe nt determina- 
tions of the stellar mas s function of galaxies (e.g. iBaldrv et al.l 
120081: lYang et al.ll2009h would lead to slightly different best-fit 
parameters. Reasonable modifications to the IMF will also yield 
slightly different best-fit parameters. However, the conclusions 
of our article are robust to these uncertainties, as they do not 
depend on the precise values of the model parameters. We also 
note that usn _ 120 km s~' is close to the value of 100 km s~' de- 
rived bv lDekel & Silkl ([1986). It implies a supernova efficiency 
of ~ 1% if there is one 10 52 erg supernova every 100 M of star 
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formation and the escape speed is y3 u c j lc , that is, the required 
wind kinetic energy is ~ 1 % of the supernova energy. 

With Li le i on = 40kms~ 1 our model cannot form galaxies be- 
low this halo circular velocity. One may argue that galaxies are 
known with lower maximum rotation velocities. Indeed, since 
our model starts at roughly the epoch of reionisation, we will be 
missing t he galaxies of lower m ass that form before reionisation 
(see, e.g jMamon et al.l2010bl) . It turns out that the precise value 
of Ureion has little effect on the conclusions of our paper, because 
of our finite galaxy mass resolution. Moreover, relative to their 
halo masses, their stellar masses are extremely low, hence they 
contribute little to the growth of more massive galaxies. 

The incorporation of ^ stl i p reduces the stellar masses of satel- 
lite galaxies at each timestep. Comparison with the case 7/ stl -i p = 
(Figure [4^, dotted line) shows that this effect is negligible at 
wJstars ~ 10 11 h~ l M Q . Tidal stripping is, however, important at 
higher masses, where the galaxy mass function is predicted to 
be shifted toward hig h masses by ~ 0.2 dex in mass compared to 
the iBell et all (120031) mass function if satellites are not stripped 
of part of their mass before they merge into the galaxies that 
populate this part of the mass function. This discrepancy may 
also be due to systematic measurement errors in the data, since 
many massive ellipticals extend outsid e the photometric aper- 
ture used to measure their luminosities. lLauer et al. I d2007h find 
errors of up to a magnitude in the SDSS luminosities of such 
gal axies and t hese are the same data that anchor the high end of 
the Hen et al.l mass function in Figure |4] We shall elaborate on 
this in the Discussion (Section l4~Tl i. 

Assuming that every halo merger immediately results in a 
galaxy merger causes the most massive galaxies to grow even 
larger (dash-dotted line in Figure |4^) but the difference between 
the dash-dotted line and the dotted line is fairly minor suggesting 
that most halo mergers do result in galaxy mergers. The differ- 
ence between the two curves in the highest mass bin disappears 
if the dash-dotted line is computed not for all halo mergers, but 
for those of bound halo-satellite systems. 

The dashed line in Figure corresponds to the extreme 
opposite assumption that whenever a new halo appears a new 
galaxy is created but there are no galaxy mergers at all. This 
assumption predicts too many low-mass galaxies and too few 
high-mass galaxies with respect to the observations. 

Both the dotted line and the dash-dotted line in Figure as- 
sume models without tidal stripping (so does the dashed line). 
In both cases the cosmic stellar mass density obtained by inte- 
grating the galaxy stellar mass function over all masses ex ceeds 
the observational value inferred from the lBell et al.l (120031) mass 
function. This suggests that stripping is necessary, since, in our 
model, changing the merging rate does not affect the cosmic stel- 
lar mass density. 

Having seen how the mass function depends on the capture 
and the probable stripping of satellites, we now examine how it 
depends on our assumptions on the quiescent growth of galax- 
ies in isolated haloes (Figure [4})). We compare the full model 
with three unrealistic simpler models, which are, however, use- 
ful for illustrative purposes. They are: i) a model in which the 
stellar mass grows proportionally to the mass of the dark mat- 
ter (m stars = /i, A^haio; dotted line); ii) a model with the term in 
fdrc but without the cutoff at high masses (dashed line); and iii) 
a model with the cutoff at high masses but not the term in v C nc 
(dash-dotted line). 

Figure 2J) illustrates the well known fact that supernova and 
reionisation feedback processes are essential to reconcile the ob- 
served galaxy mass function with the halo mass function of a 
cold dark matter Universe (compare the solid line and the dash- 



dotted line). The suppression of star formation at Mh a i > M s h oc k 
is necessary to reproduce the break in the galaxy stellar mass 
function at ~ lO n M (compare the solid line and the dashed 
line). It is the combination of these low and high mass cut- 
offs that allows the best-fit model to reproduce the normalisa- 
tion and the characteristic break of the galaxy stellar mass func- 
tion. Figure |4j} shows that the low u c j lc cut-off is important, not 
only for the low-mass end, but also for the high-mass end of the 
galaxy mass function, because the high-mass cut-off alone is not 
enough if the building blocks of giant ellipticals have been al- 
lowed to grow in an uncontrolled fashion before crossing the 
critical massQ We also find that the quenching of high mass 
galaxies by shock heating is much more effective than the reduc- 
tion of their masses by tidal stripping, as can be seen by compar- 
ing Figs. [4^ (red dotted curve) and|4j) (green dashed curve). 

Figure [5] shows the evolution of the cosmic SFR density. 
There are two ways to compute the cosmic SFR density: by ob- 
servations (or SAM simulations) of the total rate of star forma- 
tion, or by taking the time derivative of the total mass integrated 
over the galaxy stellar mass function. Since our model does not 
incorporate the concept of SFR, we adopt the latter method, 
which we show as the black lines (solid for our best-fit model). 
The agreement between our best-fit model and t he data (espe- 
cially those derived with the second method, e.g. IWilkins et al.l 
l2008l red curve) is reasonable given how simple our model is^. 

3.2. Central and satellite galaxies: the role of the 
environment 

Figure [6] provides a deeper insight of what happens inside our 
model. We have plotted all our simulated galaxies, at z = 0, in 
an m stals vs. MjjlS diagram, where M{j°j^ is the mass of the largest 
halo in which a galaxy is contained (so that if a small galaxy is 
located near a larger one, itself within a group of galaxies, the 
host halo is that of the group and not that of the larger galaxy). In 
other w ords, Figure |6] constitutes a representation of the Tully- 
Fisher (fl977t) scaling relation between stellar mass in galaxies 
and their halo properties (where the halo is quantified here by 
the mass at the virial radius instead of the maximum circular 
velocity). The symbol with which a galaxy is plotted has been 
sized and colour-coded according to the values of the fraction 
/merg of present-day mass acquired via mergers (symbol size) 
and the fraction g m \i of the mass acquired by mergers that was 
acquired in 'cold-mode' mergers (symbol colour), where 'cold- 
mode' and 'hot-mode' mergers are defined as follows. A merger 
that takes place between the two N-body snapshots s and s + 1 is 
classified as being 'cold-mode' if the mass of the host halo of the 
merger remnant at s + 1 is Mjj°^ < M s h oc k- Figure|6]only shows 
the best-fit model, corresponding to the solid line in Figure [4] 

We clearly see two galaxy populations separated by an 
empty zone: a population of central galaxies (galaxies for which 
■Mhaio = Mhak>), which follows a narrow curved band in the m stals 
- plane, roughly following the stellar masses predicted 

from the accretion channel (eq. O) and a population of satel- 
lite galaxies (Mhaio < which are scattered over a broad 
area of the diagram lying below the central galaxy relation. 



7 This result agrees with what N. Katz finds from hydrodynamic sim- 
ulations (private communication). 

8 The discrepancy in Figure [5] between the directly measured SFR 
densities (symbols with error bars) and the evolution of the SFR density 
computed from the time derivative of the integrated galaxy stellar mass 
function (red curve and shaded area) is an open observational problem. 
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Fig. 6. Galaxy stellar mass vs. host halo mass at z-0 for the set of parameters (Table[TJ that fit best the galaxy stellar mass function. 
Each point is one of our simulated galaxies. The symbol sizes and colours vary with the fraction / mer g of mass acquired in mergers, 
and the fraction g zo \& of this mass acquired by cold-mode mergers (where the sum of the masses of the haloes of the merging 
galaxies is below M s h oc k), respectively, according to the legend. The vertical green line corresponds to Mhaio = ^shock, while the 
oblique magenta line indicates m stms = f\, Mh a i - The green curves show the stellar masses of galaxies acquired through the accretion 
channel (eq. (8j with best-fit parameters from Table[TJ at z=0 (lower) and z = 3 {upper curve). The shaded cya n and gold region s 
correspond to the observations from the SDSS (Yang et al. 2009), and the halo-galaxy abundance matching of iGuo et~al 1 (I2010h . 
For the lat ter, the halo ma sses were s hifted by log 10 h to conform to our units. Also, the stellar masses were shifted by - log 10 h and 
log 10 h for lYang et al.l and lGuo et all respectively, to take into account the different units with h. We also shifted th e stellar ma sses 
by an addi tional factor of +0. 1 and +0.15 dex to pass from the Kroupa and Chabrier IMFs, respectively assumed by lYang et a"fl and 
IGuo et all to the diet Salpeter IMF assumed by Bell et al. (2003), on whose mass function we have calibrated our model. 



It is important to note that the masses in the (/-axis of Figure|6] 
are stellar masses. Neglecting the mass in cold gas is reasonable 
for mass ive galaxies, but no t at low stellar masses (log hm sVdrs < 
8.7, e.g. IZhang et ai1l2009h . where the HI mass dominates the 
stellar mass. The inclusion of gas in our toy mode l (eq. IfTOl ) 
with gas-to-star ratios given in iBaldrv et al.l d2008l and refer- 
ences therein) leads to baryonic masses of low-mass galaxies 
that are nearly twice as large as inferred from the stellar masses 
in Figure [6] 

Central galaxies with / me rg ~ 0. 1 tend to accumulate in a nar- 
row zone of the m stars - Mjj°?J diagram, which appears as a black 
curve of points, and corresponds to the present-day m^ s - Mhaio 
relation (lower green curve) given by Eq. ([8]). Its broadness is due 
to the higher value, for given M^°**, of u c j rc and therefore m stars 
at higher redshifts, as can be seen from the two green curves in 



Figure[6] Galaxies that lie well above the critical z=0 curve either 
live in haloes that have lost mass due to tidal stripping (dots and 
triangles) or are galaxies that have grown above the relation by 
mergers (filled circles). Galaxies lying below the critical z — 
curve are the centrals of hosts containing a non-negligible frac- 
tion of their mass in satellites, so that the host halo mass is larger 
than the mass of the halo directly associated with the galaxy. 

Our general m stals - Mjj"? relation predicted for the central 
galaxies of groups and clusters in Figure [6] is hig hly consis- 
tent w ith that found in the SDSS shown in Fig. 4 of lYang et al.l 
J2009I) (shaded cyan region in Fig. [6). At lower masses , it is also 
highly consistent with the relation that lGuo et al deduced 
(shaded gold region in Fig. [6]) by matching the halo and sub- 
halo mass function measured in the Millennium-II c osmological 
simulation to the galaxy mass function deduced by iLi & Whita 
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Fig. 5. The evolution of the cosmic SFR density with redshift 
in the observed Universe (points with error bars and hatched 
region) and in our model (thick black solid line; the dotted 
and dotted-dashed lines correspond to the two models where 
baryons trace dark matter and without a low circular velocity 
cutoff, respectively, see Fig. IS)). The blue sq u ares are UV data 
iGiavalisco et all (12004]); IWilson et al.l (120021): iMassarotti et al ' 
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Bunker et aH(l2004 : lOuchi et al.l (12004. The blue pentagon s 



are Hubble Ultra Deep Field estimates dThompson et al. 20061). 
The red squares are Ha, H/S and Oi l data ( Hanish et al.l 
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et al 



- .Sadler et all ._ 
iMachalski & Godlowskil l2000t lHaarsma et al.l l2000t 
Condon] |1989l). The orange pe ntagon at z 
data dGeorgakakis et all 120031). Th e green 
frared dPerez-Gonzalez et al.l 120051: iFlores et al.l 1 19991) and 
sub-mm (Bar ger et al.1 120001) data. The green hatched region 
is the far infrared (24 pm) SFR history from iLe Floc'h et al.l 
d200l. Most o f th em were taken from the compilations in 
Hopkins (2004]) and lHopkins & Beacoml d2006l) . The red curve, 
p stars = (0.014 + 0.11 z)h/[l+Sz21£^h/QJ)M Mpc- 3 , 
is the cosmic SFR derived by IWilkins et al.l d2008l) by taking 
the derivative of the cosmic stellar mass density, i.e. from the 
variation of the observed stellar mass functions with z rather 
than from star formation rate measurements. The red shaded 
area shows the margins of uncertainty around the red curve. All 
observational data have been corrected (-0.15 dex) to the diet 
Salpeter IMF assumed by Bell et al. (2003), on whose mass 
function we have calibrated our model. 



0.24 is X-ray 
squares are in- 



d2009l) . This gives us some confidence that our model of equa- 
tion (0 provides a reasonable approximation of the true stellar - 
halo mass relation. 

Figure [7^ shows how different environments contribute to 
the galaxy stellar mass function at z=0. The galaxy mass func- 
tion in haloes with 1O 12 M < hM b . d \ < 1O 123 M contains a 
strong peak at m stars ~ 10 107 h~ 1 M Q , the characteristic mass of 
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Fig. 7. The galaxy stellar mass function at z-0 split by environ- 
ment, defined by host halo mass (top) or group richness (number 
of companions with ra stars > 1O 1O M //! in the same host halo, 
bottom). The dashed curves in the top panel show the conditional 
mass function fits to the SDSS data by Yang et al. (2009). The 
same h and IMF corrections as in Fig. |6]apply. 



the central galaxy of a halo in that mass range (Fig. [7^, black 
curve). The contribution of satellites to the galaxy mass function 
is at least an order of magnitude smaller than that of the cen- 
tral galaxies. In a halo with 1O 129 M < /jM ha io < 1O 132 M , 
the relative contribution of satellite galaxies is higher but we 
still see a clear valley between central and satellite galaxies at 
fw sta rs ~ 10 104 /i _l M o , while the central galaxy peak has moved 
up to m stars ~ 1O 1U /2~'M (blue curve). The peak of the green 
curve is at a lower number density than the peak of the blue 
curve, because there are less central galaxies of haloes with 
1O 12 9 M < /iM halo < 1O 13 ' 2 M than there are central galaxies 
of haloes with 1O 12 M < hM halo < 1O 12 ' 3 M . At 1O I3 5 M < 
hMhzXo < 1O 13,8 M , the satellite population has become so nu- 
merous and the central galaxy peak has moved down so much 
that the valley is no longer visible (green curve). Thus, the more 
massive the halo, the greater is the relative importance of the 
satellite population. 

Similar bimodal (wide satellite + peaked central) galaxy stel- 
lar mass functions at give n halo mass are als o seen in the anal- 
ysis of SPSS gal axies by lYang et al.l d2008l) (their Fig. 2) and 
lYang et aD, d2009l) (their Fig. 4), with peaks and troughs occur- 
ring at very similar values of m^- drs (dashed and solid curves in 
Fig-Ek, respectively for the SDSS and our model, where we used 
the same narrow bins of halo mass as they did for better compari- 
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Fig. 8. Median fraction of the z-0 galaxy mass acquired via 
mergers rather than via gas accretion for various models: best-fit 
full model (black solid line), baryons trace dark matter (m stals = 
fb Mialo, red dotted line), no high-mass cutoff (M s h oc k — > °°, 
green short-dashed line), no low circular velocity cutoff (l>sn = 
freion = 0, blue dash-dotted line), and no tidal stripping C?7strip 
0, black long-dashed line). 

son). The lack of galaxies with stellar masses just below those of 
central galaxies (Figure|6]i thus appears to be con sistent with the 
two population (centrals + satellites) model (e.g jYang et all) . 

For observers, it is easier to define the environment of a 
galaxy in terms of the number of companions above a lumi- 
nosity limit than by its halo mass. For this reason, we show in 
Figure |7J) how the galaxy stellar mass function differs between 
rich and poor systems, where rich and poor are defined by hav- 
ing respectively less than four or at least four companions with 
wJstais > 10 l0 h~ l M o . Since we use two wide bins of richness in 
Figure|7j5, we cannot see peaks as we did in Figure^. However, 
there is an important difference between the mass functions 
of poor and rich systems: poor systems display a strong lack 
in massive galaxies relative to richer systems. In contrast, the 
low-end slopes of the mass function is the s ame for poor and 
rich haloes. In other words, if one used the ISchechterl (fl976h 
parametrisation, f(m) oc m a exp(-m/m»), the exponential cut- 
off of the mass function, m,, would be much lower for the poor 
systems, while the faint-end slope, a would be independent of 
richness. 

3.3. The importance of mergers 

Figure |8] shows the median value of the mass fraction, / merg , ac- 
quired by mergers (instead of by gas accretion) in different mass 
bins for the five models plotted in Figure^. This median hides a 
large scatter, as one can see from the distribution of f ma - g at con- 
stant m stars in Figure [6] but several statistically significant trends 
emerge, as we shall now see. 

Our standard model (black solid curve) shows that, for fi- 
nal galaxy masses above 10 h~ l M Q , typical galaxies have ac- 
quired the bulk of their mass by galaxy mergers, whereas at 
lower masses, typical galaxies grow mainly by gas accretion. 

Could the lack of mergers at low mass be a consequence 
of our limited mass resolution? One can compare our stan- 
dard model to one where galaxies follow the haloes: m stals = 
A A^haio (red dotted curve). The median fraction of mass ac- 
quired by mergers in the model where baryons trace the dark 



matter decreases with decreasing final galaxy mass for m stars 
10 10 8 /i _1 M o , which suggests possible effects of mass resolution 
at lower final galaxy masses. However, for m stars > lO ia6 /r'M , 
this decrease is much sharper in our standard model. In other 
words, as one goes from the highest galaxy masses down to 
"istai-s > 10 l0 ' 6 /i _1 M o , the median mass fraction acquired by 
mergers decreases much faster in our standard model than in 
our model where baryons trace the dark matter. We thus con- 
clude that the decreasing importance of galaxy mergers from 
'"stars = lO u -°h~ l M Q to m stals = l0 l0 - 6 h~ l M e is robust to the 
effects of our limited mass resolution. 

Note that the host haloes of galaxies with ra stars ~ 
10 105 h~ l M Q are well resolved. It is their merging histories 
that are not. Neglecting all mergers with sub-resolution haloes 
(A^haio < 1-5 x 10 9 /z _1 M Q ) does not affect the value of / merg for 
a cluster galaxy, but the effects of neglecting mergers with sub- 
resolution haloes propagate to galaxies with masses up to nearly 
the critical mass, m alt = /bM s h oc k = 1.3 x 10" lr l M Q . 

Reionisation and supernova feedback (the terms in v Q \ m in 
Eq. [8]l produce an effect analogous to that of mass resolution 
by suppressing galaxy formation in all haloes below u c i rc = 
40 km s _1 and by considerably lowering the masses of galax- 
ies in haloes just above this threshold. Therefore, mergers with 
small haloes make no or a very small contribution to the growth 
of the galaxy stellar mass, which results in a strong suppression 
of /mere at all masses (Figure [8] green dashed line vs. red dotted 
line). 

The suppression of star formation at halo masses greater than 
Mihock has a small effect on the merging histories of galaxies 
with Mhaio ~ Mshock (blue dash-dotted line vs. red dotted line in 
Figure [§), but it means that mergers provide the only mecha- 
nism for galaxies to grow above the limit mass m stals = m cr i t . 
Therefore, including the cut-off at M s h oc k increases the frac- 
tional importance of mergers in the growth of galaxies with 
Mstars > 'Went (dash-dotted line relative to dotted line), not by af- 
fecting the merger rate but by decreasing the importance of gas 
accretion. Combining the terms in u c j rc and M s h oc k in equation ^ 
causes: i) a strong suppression of the importance of mergers at 
Mstars < m cr i t and ii) a strong increase in the importance of merg- 
ers at m stals > m cr i t , both with respect to a simple model where 
the baryons follow the dark matter (dotted line). This explains 
the presence of two growth regimes, one dominated by gas ac- 
cretion, the other dominated by mergers, respectively below and 
above m stals = m cl -i t . Adding tidal stripping (solid line) has a mi- 
nor effect on the median value / melg , which may not be statisti- 
cally significant. 

Figure [6] shows that /mere increases with m stars along the re- 
lation for central galaxies. The top panel of Figure [9] splits the 
mass accretion history of galaxies into three channels: 1) gas 
accretion, 2) cold-mode mergers, and 3) hot-mode mergers (see 
Section [3~2l for a precise definition of cold/hot-mode mergers). In 
turns out that half of the galaxies with m stars > \0 n h~ l M Q have 
acquired at least half of their mass via merging, and this merging 
involves almost always a primary halo more massive than M s h oc k- 
On the other hand, most galaxies with m stars < \Q n lr l M Q have 
acquired most of their mass through gas accretion. This figure 
is virtually unchanged when we restrict ourselves to the central 
galaxies. 

The top panel of Figure|9]shows that cold-mode mergers ap- 
pear to contribute little to the growth of galaxies. In no mass 
bin is their median contribution to the final stellar mass more 
than ~ 5%. This minor channel of galaxy mass growth reaches 
its maximum relative significance at a galaxy stellar mass of 
order m stars = fi, M s h oc k- In particular, among the satellites in 
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Fig. 9. Top: Median mass fraction at z=0 acquired via gas accre- 
tion (black line), "cold mode" galaxy mergers within haloes with 
A^haio < M s hock (blue line), or "hot-mode" ones within haloes 
with M{j°^ > M s hock {red line). Error bars show the interquar- 
tile range. The green vertical dashed line shows the baryonic 
mass f7i cl -; t = ^M S hock- Middle: Same as top panel, but for galax- 
ies with /merg > 0.4 (comparable to ellipticals, see next panel). 
Bottom: Fraction of galaxies with / merg > fm for / m = 0.4, 0.5, 
and 0.6. The dark red circles show th e fracti on of ellipticals, 
as visually classified by iFukugita et alj (l2007h from SDSS im- 
ages, as a function of stellar mass, as presented bv lBernardi et alj 
(1201 Ol middle right panel of their Fig. 12). 



groups and clusters with Mjj°?J > 10 13 h~ l M Q , only 2% have ac- 
quired over half their mass by mergers (6% of the more mas- 
sive ones with m stals > 10 l0 /j _l M o and 1% of the ones with 
10 9 M o </2m stars < 1O 1O M ). 




Fig. 10. Distribution of fraction of the z-0 mass acquired by 
mergers for galaxies in different bins of stellar mass. 



Still, some low-mass galaxies do acquire most of their their 
stellar mass via mergers. Figure [6] shows nearly one hundred 
central galaxies with host halo masses below M s h oc k that have 
acquired most of their stellar mass by mergers, more precisely 
by cold-mode mergersQ Figure [TO] shows the distribution of 
the fractional contribution of mergers to the final stellar mass 
for galaxies in four stellar mass intervals. A few galaxies that 
have accreted a significant fraction of their mass via mergers 
are present even in the lowest mass bin. In other words, by 
assuming that a large mass fraction acquired by mergers im- 
plies elliptical galaxy morphology (either through single ma- 
jor mergers or throug h repeated minor mergers, as suggested by 
iBournaud et al ]2007h . our model can accommodate the presence 
of 1O 1O5 M galaxies with early type morphologies. However, 
the probability that a galaxy has an early type morphology 
(high /merg) is much higher at m stars > 10 1 1 /i 1 M than it is 
at 10 10 M o < hm staIS < 1O 1O5 M . As seen in the middle panel 
of Figure [9] cold-mode mergers thus do play an important role 
in the formation of those intermediate-mass galaxies that mainly 
grow by mergers, and which may be compared to intermediate- 
mass ellipticals. 

Finally, the bottom panel of Figure [9] shows that the frac- 
tion of galaxies where the mass fraction acquired via merg- 
ers is /merg i 0.4 matches reasonably well the observed frac- 
tion of ellipticals for galaxy masses above lQ w hr l M Q At 
IOVMq < m stars < lO'VMo, our model seems to pre- 
dict less ellipticals than observed in the SDSS. Our interpreta- 
tion (based also on the structural properties of dwarf ellipticals 
and dwarf spheroidals) is that these objects were not formed by 
mergers, but this occurs below our estimated mass resolution 
limit. We defer further discussion of this point to the conclusion. 



9 The large symbols for galaxies with high merger fractions in 
Figure [6] creates an illusion of the dominance of cold-mode mergers in 
low mass (10 < log /? m stars < 10.5) central galaxies, which is dispelled 
in Figures [9] (top panel) andllOl 

10 The drop in the observed fraction of ellipticals at very high mass 
appears is caused by an increase of SOs, but the distinction between 
ellipticals and lenticulars is difficult and subject to errors. 
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Fig. 11. Top: Fraction of galaxies with over 40% of their z=0 
mass acquired by mergers versus stellar mass in bins of halo 
mass (labelled by log hM). The errors are binomial. Bottom: 
Same versus halo mass in bins of stellar mass (labelled by 

log /!ffJ stars ). 

The top panel of Figure QT| shows the fraction of galaxies 
with over 40% of the mass acquired by mergers (as the dotted 
line of the bottom panel of Figure|9]l, split in bins of halo mass. It 
is clear the halo mass plays only a minor role in determining the 
fraction of galaxies with at least 40% of their mass acquired by 
mergers. In contrast, the bottom panel of Figure QT| shows how 
the fraction of galaxies with over 40% of the mass acquired by 
mergers varies with halo mass in bins of stellar mass. If / mere > 
0.4 is a proxy for elliptical galaxies (as can be inferred from 
the bottom panel of Figure [9j, we predict that the fraction of 
ellipticals should depend more on the stellar mass of a galaxy 
than on its global environment. 

The top panel of Figure Q~2] shows the contribution of ma- 
jor mergers to the total mass fraction at z-0 acquired via merg- 
ers for galaxies with / mer g > 0.4 (hereafter, candidate elliptical 
galaxies, following the bottom panel of Fig. |9}. Defining major 
mergers as those with stellar mass ratios between 1:1 and 1:3 
(red lines in Figure IT2T >. we find that at m stars < l0 l0 - & h~ l M e 
major mergers account for little of the mass growth of can- 
didate ellipticals occurring through mergers. But in galaxies 
with OT stals > \0 u h~ l M Q , major mergers contribute half of the 
mass growth by mergers of candidate ellipticals. If we relax 




10.5 11 11.5 

Fig. 12. Top: median fraction of the mass at z-0 acquired by 
mergers due to major mergers for galaxies with f mag > 0.4. The 
red and blue lines use stellar mass ratios between 1:1 and re- 
spectively 1:3 and 1:5 to define major mergers. The error bars 
are bootstrap estimates of the uncertainties of the mean fractions. 
Bottom: geometric mean mass ratio of mergers for galaxies with 
/merg > 0.4 at z = 0. The error bars are bootstrap estimates of 
the uncertainties of the geometric mean ratios. The red and blue 
horizontal lines respectively represent the 1:3 and 1:5 limits of 
major mergers. 



our definition of major mergers to include mergers with stel- 
lar mass ratios between 1:1 and 1:5 (blue lines in Figure \V2\. 
80% of the merger mass growth of massive candidate ellipticals 
(m > 10 11 h~ l M Q ) occurs through major mergers, while at lower 
masses (10 < log h m < 1 1), major mergers account for at least 
one-third of the mass growth via mergers. 

The bottom panel of Figure [12] shows the geometric mean 
of the mass ratios of the merging galaxies involved in the mass 
growth of candidate ellipticals. Again one sees that the merger 
mass growth of candidate ellipticals is mainly through minor 
mergers at low mass (< 10 105 h~ l M Q ) and major mergers at high 
mass (> 10" h~ l M Q ). The importance of major mergers peaks 
around f,z stars = 10 u - 3 /z _1 M o . 



4. Discussion and conclusion 

Our goal is to understand how galaxies grow in mass, i.e. what 
are the respective roles of gas accretion and galaxy mergers, as 
well as that of feedback mechanisms in the growth of galaxies. 
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We have developed a toy model of galaxy formation, using 
a very simple hybrid SAM/HOD approach to parameterise, as 
a function of halo mass and redshift, the stellar mass present 
in galaxies after gas accretion and quenching of gas infall by 
virial shocks and AGN and of star formation by reionisation and 
supernovae. We also include galaxy mergers, as measured from 
our high-resolution cosmological N-body simulation. The stellar 
mass acquired via gas accretion is modelled with a simple equa- 
tion involving only three parameters, while we added a fourth 
parameter to describe the stellar mass loss caused by tidal strip- 
ping. We fined-tuned our four parameters to the observed z-Q 
stellar mass function of galaxies. 

In contrast with SAMs, our toy model does not incorpo- 
rate structural or morphological properties of galaxies, has no 
time derivatives built-in (for the transfer of mass between galaxy 
components), and does not incorporate luminosity, hence our 
galaxies have no colours. 

Our choice to model the stellar mass as a unique function 
of halo mass and redshift may seem oversimplified in view of 
the complexity of galaxy formation, as exemplified by the very 
complex semi-analytical codes of galaxy formation. However, at 
log/iMhaio > 12.15 (where it can be measured) the SDSS ob- 
servations also appear to show a very smal l scatter of m stais vs. 
Mh- d \o for central galaxies dYang et al.l l2009 and the cyan region 
of our Fig. |6j». Moreover, both semi-analytical and hydrodynam- 
ical models of galaxy formation indicate that central galaxies 
in halos show a ver y small scatter betwee n the stellar and halo 
masses (see Fig. 2 of lCattaneo et al.ll2.Q07b , at z — and at z — 3. 

The main worry with our model is that our prescription for 
wJstars of equation © may not be correct at z > 0. However, our 
m stars -Mhaio relation at z — 5 (Fig.pt is not dramatically different 
from that of iBehroozi et al.l d2010l) . Moreover, we have reasons 
to believe that the cutoffs at low and high mass should be of the 
right order, so the worry is on the maximum efficiency (e.g. the 
maximum height of the top [red] curve of Fig. [2j. Furthermore, 
our approximate match of the evolution of the cosmic star for- 
mation rate density (Fig. [5J suggests that, at z > 0, our model of 
equation © is of the right order of magnitude. 

While we do not claim to have explored all possible parame- 
ter combinations, any model that suppresses star formation even 
further will worsen our agreement with the normalisation of the 
cosmic SFR density (Figure [3}. We also note that the tendency 
to be at the lower limit of the observational range for the cosmic 
SFR density is not specific to our model, since it is also found in 
state-of-the-art semianalytic models that are based on the sam e 
fundamental assumptions (see Figure 8 in lCattaneo et al.ll2006l) . 

4.1. Tidal stripping 

Without tidal stripping, we find that, even after correcting for 
overmerging, we still have an e xcess of massive galaxies com- 
pared to the iBell et all (120031) mass function (dotted line in 
Figure |4^). This remains true for other o bserved determina- 
tions of the gala xy stellar mass function ( iBaldrv et al.l l2008t 
lYang et a"Dl2009l) . One could argue that, although the semian- 
alytic model that we use for the merger tim e is a good matc h 
to what is seen in hydrodynamic simulations dJiang et al.ll2.Q08l) , 
we may still overestimate the galaxy merger rate, which is still 
somewhat uncertain. However, it is difficult to lower the merger 
rate without overproducing low-mass galaxies (see the dashed 
line in Figure |4^). One could imagine compensating for this ef- 
fect by increasing feedback at low masses, but we could not find 
a plausible parameter set that does this and fits the galaxy mass 
function as well as the full model in Figure [4] does. 



We believe that tidal stripping and observationa l errors in 



the SDSS magnitudes that anchor the high end of the 



d200l mass function (which according to lLauer et al 



Bell et all 



2007 can 



run up to a magnitude) can easily reconcile our model with the 
high end of the galaxy stellar mass function. The fundamen- 
tal reason why the magnitudes of bright ellipticals are difficult 
to measure is that their shallow light profiles do not converge. 
So it is hard to estimate the total light, let alone distinguish 
galaxy light fro m intracluster light ( see the discussion in the 
Appendix A of ICattaneo et al.l |2008h . Since the extended en- 
velopes of massive ellipticals are probably the products of inter- 
actions, the distinction between the full model and the full model 
without stripping may be one more of name than of substance. 
If we decide that the envelopes belong to the galaxies, then we 
should conclude that the full model without stripping is correct 
and that the problem lies with the data, which are likely to under- 
estimate the real masses. If instead we decide that the envelopes 
are intracluster light, then we should conclude that the model 
that includes tidal stripping is the most physical one. 

Our best-fit stripped fraction per orbit, 77 stl -; p = 0.4, matches 
well the value determined in a recent dissipationless simula- 
tion of a dwarf i rregular galaxy orbiting th e Milky Way with 
a fixed potential dKlimentowski et aljfcOOgh Fl For T] sbip = 0.4, 
the stripped mass in haloes with M ha i ~ 10 13 - 10 14 /i _1 M o is in 
the range of ~ 0.05 - 0.1/bMh a i o - This v alue is consistent with a 
numbe r of independent measurements: ICastro-Rodriguez et al.l 
(12009b find 7 % of diffuse intr a cluster light (ICL) in the Virgo 
cluster, whil e Feldmeie r et al ] d2004h find 15% in Virgo; but 
Lin &Mohrl (12004 find as much as 50% in clusters, suggesting 
that ?7 stl jp might be even higher than we found. 

4.2. The role of feedback 

In agreement with previous studies, we find that in all haloes, 
but more so in low (M ha i <K M shock ) and high (M ha i > M shock ) 
mass haloes, m stars /Mh a i must be much smaller than the univer- 
sal cosmic baryon fraction if the halo mass function predicted 
by the CDM model is to be reconciled with the observed galaxy 
mass function. 

The mechanism that suppresses star formation in haloes with 
A^haio > Mshock is clear. The post-shock gas is hot enough to 
maintain a stable shock, while in lower mass halos, the post- 
shock gas radiates efficiently and can no longer provide the ther- 
mal pressure to support the shock and collapses instead. The 
question is why can't the gas around massive haloes cool down 
again. Growing evidence suggests that this is due to heating 
from the central AGN (see ICattaneo et al.ll20 09; and references 
therein). 

Reionisation and supernova feedback are widely invoked to 
explain the very small OT stars /Mh a i ratios of galaxies in low mass 
haloes (e.g.. lMcGaugh et"a"fll2010l) . The trouble is that attempts 
to simulate these processes have difficulty to produce outflows as 
large as those that are required by this study because the energ 
that is available is used i nefficiently dDubois & Tevssierll200 
ICeverino & Klvpinll2009l) . However, unless the CDM model is 
wrong, whatever other process may be relevant must behave 
analogously to the reionisation and supernova feedback model 



11 In a sequel with a live Milky Way potential, [Lokas et al.l 1201 Oh find 
an even stronger decrease in stellar mass, that amounts to 0.65 per or- 
bit among the particles within 1 kpc. Note, however, that the secondary 
galaxy is 4 times more massive an d with half the apocenter as the re- 
spective values in the simulations of Klimentow ski et alj d2009h . 
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describe d in the article, i f it is to fit the observed stellar mass 
function (iBell et al.ll2003l) that we have used. 

The fact that m stars /Mh a io drops at halo masses below 
and above M s h oc k introduces a characteristic scale for galaxy 
masses that mak es them deviate from the dark matter's self- 
similar evolution jMarinoni & Hudsonll2002l : iBaldrv et aLll2Q08t 
iLi & Whitel 12009-1) . If the evolution of baryons followed that of 
the dark matter, then galaxies of all masses would have similar 
merging histories (to the extent that dark matt er halos evolve 
self-s imilarly, which is true to first order, e.g. Ivan den Boschl 
120021) . Instead, it is clear that the baryonic physics, whatever 
they are, break the dark matter's scale-invariant behaviour by 
suppressing star formation in low and high-mass haloes. It is 
only within haloes in a narrow mass-range around Mhaio ~ 
M s hock ~ lO 12 /i _I M that baryons can form stars efficiently (the 
galaxy formation zone in Fig. |3). 

This conclusion is similar to that of iBouche et al] (1201 Oh . 
who developed a toy model similar to ours in that it also incor- 
porates minimum and maximum halo masses for efficient galaxy 
formation, although their emphasis is different. They argue that 
a fixed halo mass floor of 1O H M for galaxy formation appears 
to be required to reproduce the observations of specific star for- 
mation rate as a function of stellar mass, the evolution of the 
cosmic star formation rate and the Tully-Fisher relation. Our 
model is more focussed on the accurate treatment of the role 
of galaxy mergers in the growth of galaxies, than on the model- 
ing of star formation rates. However, our model also allows for 
lower galaxy masses at high redshift, which may help explain 
the lowest mass dwarf galaxies. 

4.3. Mergers vs. gas accretion 

The shutdown of gas accretion above the critical mass M s h oc k ~ 
10 12 M o is the reason why, in our model, mergers are the only 
opportunity for galaxies to grow above the limiting mass m cr i t = 
A Mihock- Therefore, galaxies with m ga i » m cr i t must have ac- 
quired a significant fraction of their mass via mergers. This con- 
clusion appears quite robust because without this sh utdown mas- 
sive galaxies would be too frequent and too b lue dBower et al.l 
l2006btlCatta"neo et al.ll2006l ICroton et al.ll2006h . 

Where does the dearth of mergers at galaxy masses lower 
than lO'^-'Mo come from? It cannot be a resolution effect 
(though finite resolution makes it more pronounced) because 
even at m stars > 10 10 - 6 li~ l M Q the mass fraction acquired via 
mergers rises faster with galaxy mass in our best-fit model than 
when we force the baryons to trace the dark matter. There is, 
instead, a simple explanation why the importance of mergers de- 
creases at low masses. The steeper slope of stellar versus halo 
mass at low mass (Fig. [6} naturally reduces the relative abun- 
dance of major partners relative to minor ones, and therefore 
the ma ss contribution that any merger can make. IHopkins et al] 
J20T0I) have made a similar kind of argument. 

The present work reinforces the new picture in which most 
galaxies grow through gas accretion and form stars along the 
Blue Cloud until the critical halo mass for shock and AGN 
heating is reached. Then, gas accretion is shut off, star forma- 
tion is quenched and they migrate to the Red Sequence, where 
gas-poor mergers represent their only opportunity for growth 



(Bower et al. 2006b; Cattaneo et al.l l2006t ICroton et al 



iFaber et alJl2007LICattaneo et al.ll2008l)T 



We refer to mergers that take place in host haloes above 
A^shock as hot-mode mergers, and to mergers that take place in 
host haloes below M s h oc k as cold-mode mergers. The former 
usually contain little cold gas, while the latter involve galaxies 



containing large masses of cold gas supplied by cold accretion. 
It is thus tempting to identify hot-mode and cold-mode merg- 
ers with gas-poor (dissipationless) and gas-rich mergers, respec- 
tively. The latter are mainly confined to galaxies below the crit- 
ical stellar mass _/bM s hock (Fig. [6j and are most frequent just 
below f D M s hock (top panel of Fig. |9j. This points to a char- 
acteristic mass of ~ IP 105 /; -1 M for u ltraluminous infrared 
galaxies (ULIRGs , e.g. lColina et al.ll2001l) and quasar hosts (e.g. 
iBonoli et al.ll2009t) . of which gas-rich mergers are the supposed 
triggering mechanism (top panel of Figure|9]l- 

Our model, therefore, predicts three regimes of galaxy 
formation (see top panel of Figure [5J: Galaxies with 
"ismrs ~ 10 10 M Q were effectively built up by gas accretion only. 
At 1O 1O M i m stars t 10 11 M gas accretion remains the dom- 
inant mechanism, but we also see a population that was built 
by gas-rich mergers. Galaxies with m stars i 1O H M were mainly 
built by gas-poor mergers, in agreeme nt with the analysis of 
SDSS galaxies by iBernardi et al.l (1201 ll) . who explain the high- 
mass trends of color, elongation and color gradient of SDSS el- 
lipticals with gas-poor mergers, by invoking a transition at the 
same stellar mass. Galaxies built by gas-rich mergers are rare at 
all masses. However, as expected from our model (middle panel 
of Fig. [5), cold-mode mergers dominate the growth of those 
galaxies with TO stars < 10 n O that are mainly built by mergers. 

Runnin g a semi-analytical mo del on top of the Millennium 
Simulation, iDe Lucia et al. 1 (120061) find that the 'effective' num- 
ber of progenitors (a proxy for the number of major progenitors) 
of galaxies increases sharply beyond 1O H M , which appears to 
be consistent with our conclusion that high-mass galaxies are 
built by major mergers. 

Analyzing the same Millennium Simulation (which has 1 1 
times worse mass resolution than our simulation). iGuo & White! 
(HH) measured the rate of growth of galaxies in the channels 
of gas accretion (which they denote "star formation"), and ma- 
jor and minor mergers. They also found that while major merg- 
ers dominate the growth of the most massive present-day galax- 
ies, mergers contribute negligibly to the growth of low present- 
mass galaxies. Moreover, they noticed, as we did (our Fig. [8} 
that the importance of mergers increases with final galaxy mass 
in a stronger way t han do mergers of the dark matter haloes. 

IHopkins et all (l2010l) followed an alternative approach to 
analyse the role of mergers in the growth of galaxies, in par- 
ticular the relative importance of major and minor mergers. 
They used the halo merger rate function that iFakhouri & Mai 
(120091) derived from N-body simulations together with abun- 
dance matching to associate stellar masses to the merging haloes. 
The delay between halo mergers and galaxy mergers wa s calcu - 
lated with the fitting for mula by Bovlan-Kolchin et"ai] d2008l) . 
which is analogous to the lJiang et al.l ( 2008 ) formula that we use 
here. Their results for the relative importance of major mergers 
are, to first order, similar to ours: 1) major mergers dominate the 
formation and assembly of m stars ~ m cr i t ellipticals, but the con- 
tribution of minor mergers is non-negligible ~ 30%; 2) the for- 
mation of elliptical gala xies becomes do minated by minor merg- 
ers at lower masses; But IHopkins et al.| find that the relative im- 
portance of major mergers rapidly decreases again at high galaxy 
masses m stars > 1O U M , while we find a small and not statisti- 
cally significant drop in the importance of major mergers at very 
high galaxy masses (log /1 m stals > 11.5). When we consider the 
geometric mean mass ratio as a function of stellar mass, which 
is more comparable to the median mass ratio of bulges consid- 
ered by IHopkins et al.l we obtain (bott om panel of Figure [T2l a 
globally similar trend as lHopkins et al.l although again with two 
differences: 1) their peak for for major mergers occurs at 3 times 
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lower galaxy mass (after correcting for h); 2) the drop in merger 
mass ratio is much more pronounced in their model than in ours. 

It is difficult to say which result i s more correct in quantita- 
tive detail. The approach followed bv lHopkins et al.l is certainly 
more accurate in its capacity to predict the galaxy merger rate at 
a given time, since the halo occupation distribution is directly de- 
rived from abundance matching. However, our approach follows 
more consistently the merger history across the Hubble time of 
eac h halo inside the N- body simulation. 

iMaller et alj d2006l) analysed the problem directly with hy- 
drodynamic cosmological simulations. They concluded that 
mergers are more important at high galaxy masses. However, 
their simulations lack effective feedback and substantially over- 
estimate the observed galaxy mass function. Therefore, they 
overestimated the role of gas accretion. Nevertheless, they deter- 
mined that the average number of both minor and major mergers 
increases with ra stals . They found that their high mass galaxies 
have typically undergone one major (1:1 to 1:4) merger through 
their lifetime. We note that the threshold they use to define high- 
mass galaxies is low (log/ira stars = 10.6), so we expect that the 
dominance of mergers in their high mass bin is diluted by the im- 
portance of accretion near this mass threshold. Not only do we 
find that mergers are common among the most massive galaxies. 
We also find that the masses of these galaxies were mainly built 
by mergers. 

iNaab et aTl d2007l) also used hydrodynamic simulations to 
address the i mport ance of mergers. However, in contrast with 
IMaller et al.l d2006l) . who simulated a cosmic volume at low res- 
olution, they used very high resolution to zoom into the forma- 
tion of three individual ellipticals. They found that their galaxies 
have grown by mergers from z=l to z=0 by 25% in mass on 
average. They, therefore, concluded that intermediate-mass el- 
liptical galaxies were not built by major mergers. We note, how- 
ever, that their three ellipticals accreted over half of their mass 
by mergers since z-5. Moreover, at z=0, these galaxies live in 
haloes with virial masses in the range of 1.6 - 2.3 x 10 12 M o . 
For central galaxies, this halo mass range corresponds to m stals 
lO llo h~ l M , which is the precise galaxy mass for the transi- 
tion in mass acquisition 1) by gas accretion and by mergers 
(Figure |9), and 2) by minor and major mergers (Fig. \T% . We 
therefore find the agreement between our results and those by 
iNaabet al.lri2007l) quite encouraging. 

The increasing importance of mergers on one hand and the 
importance of major mergers on the other are confirmed through 
two analytical developments that are deferred to a forthcoming 
paper. 

Finally, we remark that the delay of galaxy mergers with re- 
spect to halo mergers, w hich we compute w ith the the dynami- 
cal friction formula in the lJiang et all (120081) version, turns out to 
play a very modest role in our conclusions: the galaxies that have 
not yet merged with the central galaxies of their direct hosts, 
even though their galactic haloes have merged with their imme- 
diate host haloes (or have become unresolved because of tidal 
stripping of most of their mass), account for i 1 % of the total 
galaxy population and ~ 10% of the satellite galaxy population 
at ra stars > l0 9 h~ l M Q (satellites make up ~ 10 - 15% of all 
galaxies). 



4.4. The dichotomy between central and satellite galaxies 

Our very simple toy model of galaxy formation shows very dis- 
tinct properties for central and satellite galaxies (see Fig. |6|. 
Central galaxy mass increases slightly with halo mass, while 
the galaxy mass function, when measured within a narrow range 



of halo mass, clearly separates the central population and the 
satellite population, with a sharp peak at high-masses that corre- 
sponds to central galaxies (top panel of Fig.[7]i. This feature was 
first predicted bv iBenson et al.l (120031) andlzheng et al.l (120051) 
and first clearly observed bv lYang et al.l (120091) in the SDSS sur- 
vey, who also modelled it with the conditional stellar mass func- 
tion fo rmalism. As we wrote these lines, we learnt that lLiu et al.l 
were also able to reproduce this behaviour with a semi- 
analytical model. 

Over three-quarter£3 of the galaxies in clusters, the huge 
majority o f which are satell ites, are observed to lie on the Red 
Sequence dYang et al.ll2008l), which is visible to very low lumi- 
nosities (e.g. iBoue et al.l2008b and is mainly composed of galax- 
ies with early-type morphologies, the great majority of which are 
classified as dwarf ellipticals (dEs). 

Since mergers are found to be unimportant for the growth 
of most satellites (Fig. [6j, including those in groups and clus- 
ters (Sect. 13.3b . one is led to the conclusion that dEs must also 
acquire their mass by gas accretion. Therefore, the morphologi- 
cal properties and the shutdown of star formation in dEs (and by 
extension in dwarf spheroidals) must be due to processes unre- 
lated to mergers and not included in our model (which does not 
describe t he conversion of gas into stars), s uch as ram-pressure 



stripping dGunn &Gottl[T97l . starvation dLarson et al 



1980) 



1998 



and harassment by repea ted fast encounters dMoore et al.1 
iMastropietro et al.ll2005b . 

Therefore, we conclude that the mass-growth of ellipti- 
cal galaxies is a function of their present-day mass: the most 
massive ones (m stars > lO 11 /^ 1 M ) have mainly grown by 
gas-poor mergers, intermediate-mass ellipticals (IO^/z^Mq < 
wistars < 10 11 h~ M Q ) by gas-rich mergers, and the low-mass ones 
('"stars < 10 I0 /i _1 M Q ) by gas accretion and later transformed by 
non-merging processes such as ram pressure stripping or galaxy 
harassment. This is similar to the picture for all galaxies, ex- 
cept that intermediate-mass ellipticals grow by gas-rich mergers 
while intermediate-mass spirals are built by gas accretion. 



4.5. Caveats 

An important caveat to our analysis is that our model does not 
describe star formation rates, as we are only interested in the 
total mass of a galaxy and the fraction of this mass that comes 
from mergers. This paper is not concerned with the fraction of 
the mass that is in gas and the fraction of the mass that is in 
stars. The SFRs that we show (Figure [5]) are simply mass ac- 
cretion rates that are estimated assuming that all the accreted 
gas is instantaneously turned into stars. Therefore, it is not sur- 
prising that our model anticipates the peak of the star forma- 
tion rate in comparison with observations (Figure |5J. In prac- 
tice, at redshift z ~ 0, where we compare our model with the 
data (Figures and [6}, the gas is usually a small fraction of the 
galaxy mass for inter mediate and high-mass galaxies (see, e.g., 
iMcGaugh et al.ll2010l and references therein). Therefore, the er- 
ror that we make by identifying the stellar mass with the total 



12 Th e 22% of non-Red Sequence cluster galaxies found by Ya ng et al .1 
(2008) matches perfectly the fraction of interlopers projected along the 
virial sphere, i.e. the average fraction of particles in the virial cone that 
are outside the virial sphere dMamon et al.ll20103) . Hence, the fraction 
of non-Red Sequence galaxies within the virial sphere is probably much 
lower than one-quarter. Indeed, the observed increase of the fraction of 
recent star for ming galaxies (RSBG s) with projected radius has been 
deprojected by Mah aian et alj feOlOh to yield a fraction of 13 ± 1% of 
RSBGs within the virial spheres in comparison with 18 ± 1% within the 
virial cones. 
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galaxy mass (stars plus cold gas) is at most ~ 10 - 20%, which 
is compatible with the level of accuracy that we expect from 
our model. We note that even observationally, the cosmic SFR 
density derived from dust-corrected measured SFRs deviates at 
z > 1 from the time derivative of the cosmic stellar mass density 
(IWilkins et alJl200§: Figure[5]l. The origin of this discrepancy is 
an open problem. Our model matches the observed time deriva- 
tive of the comic stellar mass density much better than it matches 
the measured star formation rate density. 

Much higher mass resolution on our side would certainly be 
desirable to properly follow substructures until they have merged 
with their hosts. However, this problem has been partly handled 
with the incorporation of delayed merging on a dynamical fric- 
tion timescale (Sect. 12.2b . Moreover, we are confident that our 
finite mass resolution is sufficient to indicate the strong drop in 
the importance of mergers at masses below IO'Vi^'Mq (Fig. [8]). 
Also, our low Veto cut-off implies that there is not much advan- 
tage at resolving haloes with masses that are much lower than 
the mass at which ii c ; rc = u re ion, because there is no formation of 
galaxies in those haloes after the epoch of reionisation. 

Of course, the details of our toy galaxy formation model 
are most probably oversimplified. The trend seen in the relation 
of galaxy mass to halo mass (Fig. [6]) might be incorrect at low 
masses (althoug h there is good agreement with the relation that 
iGuo et al"1l2010l derived from abundance matching, see Fig. [6]). 
Nevertheless, none of our conclusions seem to depend on the 
accuracy of the model at low masses. 

4.6. Open questions 

Mergers are a reality and we have tried to present a careful 
and robust estimate of their importance. Our work suggests that 
present-day giant ellipticals can only be built by gas-poor merg- 
ers, as gas accretion must be quenched at high halo mass to 
avoid over-massive galaxies. Moreover, intermediate-mass ellip- 
ticals have probably grown by gas-rich mergers, while mergers 
appear unimportant for dwarf ellipticals. These conclusions are 
consistent with the structural properties of ellipticals: the sur- 
face brightness profiles of massive ellipticals tend to display in- 
ner cores, while those o f less massive ellipticals tend to be cuspy 
(iKormendv et al.ll2009l) . 

Still many questions remain. Can the well-defined relation 
between the masses of ce ntral supermassiye blac k holes and the 
spheroids of their hosts dMagorrian et al.l["l998l) be understood 
with both gas-poor mergers at the high mass end and gas-rich 
mergers at intermediate masses? Can this relation be extended 
from ellipticals to spiral galaxies? How do mergers fit in the 
growth of galaxies viewed at much higher redshift? Are other 
mechanisms resp onsible for the gr owth of proto-ellipticals at 
high redshift (e.g. lDekel et"aT1l2009h ? 
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